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ABSTRACT 

The  approximate  analysis  of  three-dimensional  incompressible 
turbulent  boundary  layer  downstream  of  vortex  generators  is 
presented.  Extensive  numerical  computations  are  carried  out  to 
assess  the  effectiveness  of  single-row,  counter-rotating  vane-type 
vortex  generators  to  alleviate  flow  separation.  Flow  separation 
lines  downstream  of  the  vortex  generators  on  a  thick  airfoil  are 
determined  in  terms  of  size,  location,  and  arrangement  of  the 
vortex  generators.  These  lines  are  compared  with  the  separation 
line  without  the  vortex  generators.  High  efficiency  is  obtained 
with  the  moderately  slender  rectangular  blade  of  the  generator. 
The  results  indicate  that  separation  is  alleviated  more  effectively 
in  the  region  closer  to  the  symmetry  axis  of  the  generator  than  in 
the  outer  region  of  the  symmetry  axis.  No  optimum  conditions 
for  the  alleviation  of  flow  separation  are  established  in  this 
investigation,  and  no  comparisons  are  made  with  other  analytical 
results  and  experimental  data. 
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INTRODUCTION 

Structural  considerations  require  the  root  cross-sections  of  large  aspect  ratio 
wings  to  be  thick.  Flow  over  such  a  thick  airfoil  section  is  characterized  by  a  rapid 
pressure  rise  toward  the  trailing  edge,  and  this  large  adverse  pressure  gradient  causes 
the  boundary  layer  to  separate  from  the  surface  of  the  wing.  The  resulting  stall  is 
generally,  undesirable,  and  therefore  it  is  necessary  to  alleviate  the  flow  separation. 
Among  the  various  methods  of  suppressing  flow  separation  is  the  use  of  vortex 
generators,  which  has  been  found  to  be  relatively  effective  and  simple. 

Many  types  of  vortex  generators  have  been  experimentally  investigated  from  time 
to  time.  Schubauer  and  Spangenberg1  presented  10  different  shapes  in  their  1960 
survey.  In  the  present  analyses,  the  flow  field  of  a  geometrically  simple  vane-type 
vortex  generator  is  investigated. 


Since  the  vortex  generator  flow  field  on  a  wing  is  three-dimensional  and 
turbulent,  the  three-dimensional  incompressible  turbulent  boundary-layer  equations 
are  analyzed  by  extensive  numerical  computation.  Several  parameters  are  investigated, 
including  the  length  of  a  vane,  the  width  between  opposite  vanes,  the  chord  and 
height  of  a  vane,  and  the  locations  of  the  vortex  generators.  The  separation  lines 
downstream  of  the  vortex  generators  are  determined  to  assess  the  extent  of  the  flow 
separation  alleviated  by  the  vortex  generators. 

A  survey  of  the  literature  has  shown  that  investigations  of  the  vortex  generator 
flow  field  are  limited.  The  following  researchers  are  among  those  who  have 
conducted  investigations  of  vortex  generator  flow  fields: 

•  Brown  (1947)2  investigated  the  induced  velocity  field  of  a  vortex  generator  and 
prevented  separation  in  the  inlet  diffuser  of  an  aircraft  engine  by  using  a 
vortex  generator. 

•  Jones  (1957)3  investigated  the  vortex  path  of  a  wing-type  vortex  generator  and 
evaluated  the  vortex  strength  of  the  generator  using  potential  flow  methods. 

•  Schubauer  and  Spangenberg  (I960)1  measured  the  boundary-layer 
characteristics  of  various  vortex  generator  flow  fields  showing  the  mixing  of 
the  high  velocity  outer  layer  with  the  low  velocity  inner  layer. 

•  Pearcey  (1961)4  and  Gartling  (1970)5  investigated  the  effects  of  vortex 
generators  on  flow  separation  alleviation  at  transonic  speeds. 

•  Gartling  (1970)s  also  studied  the  vortex-generator  flow  at  supersonic  speeds. 

•  Senoo  and  Nishi  (1974)6  carried  out  experiments  for  conical  diffuser  flow 
separation  alleviation  at  divergence  angles  of  8,  12,  16,  20  and  30  deg  for  an 
area  ratio  of  4  using  a  variety  of  blade  arrangements,  inlet  boundary-layer 
thicknesses  and  locations  of  vortex  generators.  The  results  showed  that  the 
vortex  generator  prevents  flow  separation  up  to  a  16-deg  divergence  angle  and 
that  the  pressure  recovery  is  approximately  equal  to  that  of  the  best 
conventional  conical  diffuser.  At  blade  setting  angles  of  14-16  deg,  the  best 
pressure  recovery  is  achieved  and  larger  aspect  ratio  blades  were  the  most 
efficient.  At  a  nominal  diffuser  divergence  angle  of  8  deg,  pressure  recovery  is 
not  affected  by  an  increase  in  vortex  generator  number;  at  16  deg,  pressure 
recovery  remains  the  same  in  spite  of  an  increase  in  the  vortex  generator 
number  beyond  10. 

•  Rousseau  (1984)*  carried  out  a  limited  experimental  investigation  of  vortex 
generators  attached  to  a  thick  airfoil. 

•  Nickerson  (1986)7  studied  vortex  generator  flow  at  low  Reynolds  numbers. 

•  Greene  (1986)8  modeled  the  decay  of  a  wing  tip  vortex  in  the  atmosphere  in  a 
wake  downstream  of  a  wing  in  the  atmosphere  and  found  that  typical  vortex 
wake  decay  is  caused  by  flow  instability  (vortex  pair  linking)  except  under 
strong  stratified  conditions. 

•  Lee  (1986)9  analyzed  two-dimensional  vortex  generator  flow  and  showed  that 
flow  separation  can  be  prevented  on  a  thick  airfoil. 
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FLOW  FIELD  AFFECTED  BY  VORTEX  GENERATORS 

Examples  of  various  vane-type  vortex  generators  such  as  co-rotating,  biplane, 
and  counter-rotating  vortex  generato:  .  are  shown  in  Fig.  1.  Of  these  vortex 
generators,  the  counter-rotating  vane-type  generators  are  considered  in  the  present 
analysis.  The  specific  arrangement  and  nomenclature  are  shown  in  Fig.  2. 

Figure  3  shows  the  line  vortices  of  strength  r  in  the  y-z  plane.  These  line  vortices 
correspond  to  the  tip  vortices  originating  at  the  tip  of  the  rectangular  plate  vortex 
generators  which  are  at  angle  of  attack  to  the  upstream  flow.  Since  the  incident  angle 
considered  in  the  present  report  is  less  than  35  deg,  the  tip  vortices  are  steady,  and 


the  vortex  lines  follow  the  streamlines  (see  Lugt,10  p.  81).  The  counter-rotating  tip 


vortex  lines  of  equal  vortex  strength  are  straight  and  parallel  to  each  other  (Lugt,10 
p.  38).  Jones  (1957)3  and  Widnall  (1975)11  also  indicate  that  the  tip  vortex  lines  are 
straight  and  parallel  to  each  other  in  the  x-z  plane  and  prevail  a  considerable  distance 
downstream.  Jones  ( 1957)3  studied  the  vortex  path  of  the  counter-rotating  vortex 
generator  and  indicated  that  the  vortex  first  approaches  the  surface  and  then  moves 
away. 

To  achieve  high  performance,  Schubauer  and  Spangenberg  (I960)1  suggested  a 
blade  height  of  hv  =  1.2  6.  Thus,  the  height  of  the  vortex  line  over  the  surface  is 
assumed  to  be  that  of  the  boundary-layer  outer  edge,  i.e.,  h  ~  6  in  this  analysis,  so 
that  a  more  effective  vortex  flow  to  alleviate  the  flow  separation  may  be  expected. 
However,  in  order  to  correlate  the  geometric  and  position  parameters  affected  by 
strong  mixing,  a  larger  value  of  h  is  also  used  for  the  computation.  The  origin  of 
such  a  vortex  line  is  assumed  to  be  the  midpoint  of  chord  cv.  Since  a  three- 
dimensional  vortex  generator  flow  field  is  complex  and  deals  with  turbulent  flow,  no 
exact  analytical  solutions  are  expected.  Hence,  for  convenience,  it  is  assumed  that  the 
direction  of  the  tip  vortex  line  downstream  of  the  generator  is  the  same  as  the 
undistrubed  flow  direction;  the  neighboring  vortex  does  not  influence  the  vortex 
path,  and  the  vortex  does  not  burst.  Comparative  studies  of  the  extent  of  the  flow 
separation  allieviation  in  terms  of  D/d,  D/h,  cv/h,  location  of  generators,  and 
without  a  vortex  generator  are  carried  out  so  that  the  optimum  conditions  of 
separation  alleviation  may  be  determined.  The  coordinate  systems  of  the  vortex 
generator  and  the  wing  are  shown  in  Figs.  3  and  4. 

The  induced  velocity  at  an  arbitrary  point  x,y,z  due  to  tip  vortex  line  A  is  given 
by  Biot-Savart  law  as 

(v4)  =  —j==J=====  [i  +  - . -  *v  -  1  (1) 

A  4nN/  (y  -  h)2  +  (z  -  d/2)2  |_  V  (x  -  xv)2  +  (y  -  h)2  +  (z  -  d/2)2  J 


where  xv  =  vortex  origin  location, 

h  =  fip  vortex  core  height  as  well  as  tip  vortex  stream  line  height, 
d  =  distance  between  two  neighboring  counter-rotating  vortex  lines, 

T  =  vortex  strength  and  its  clockwise  direction  (positive  when  viewed  from 
front  to  back). 

The  induced  velocity  in  the  z-direction  is 
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x  -  xv 


(x  -  xv)2  +  (y  -  h)J  +  (z  -  d/2)2 


/,v,  / 
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Then  the  induced  velocity  in  the  z-dircction  due  to  tip  vortex  line  A.B.C.D 
becomes 

w _ _  h~y  [.  ^ _ x~xv _ ] 
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Thus,  the  induced  velocity  in  the  z-direction  due  to  the  vortices  shown  in  Fig.  3  is 


(y  -  h)2  +  [z  -  (kD  +  d/2)]2  \  V  (x  -  xv)^  +  (y-h)^+[z- (kD  +  d/2)] 


(y  -  h)2  +  [z  -  (kD  -  d/2)]2  {  V(x-  xv)2  +  (y  -  h)^  +  [z  -  (kD  -  d/2)] 


(y  +  h)2  +  [z  -  (kD  -  d/2)]2  1  V(x-xv)^  +  (y  +  h)^  +  [z  -  (kD  -  d/2)] 


(y  +  h)2  +  [z  -  (kD  +  d/2)]2 
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The  complexity  of  the  generator  flow  field  and  the  low  aspect  ratio  of  the  blade 
make  accurate  evaluation  of  T  a  formidable  task.  Thus  T  is  evaluated  using  the  Jones 
(1957)3  approximate  solution  of  strength  T  based  on  the  lifting  line  theory;  see 
Appendix  A.  For  a  given  x  and  z, 
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where  ARC  =  cv/h. 
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may  be  evaluated  using  Eq.  4. 


BOUNDARY-LAYER  FLOW  ANALYSIS 

If  the  boundary  layer  is  thin  and  the  Reynolds  number  is  large,  the  steady  three- 
dimensional  boundary-layer  equations  for  turbulent  incompressible  flow  may  be  used. 
With  the  coordinate  system  shown  in  Fig.  4  these  equations,  taken  from  Nash  and 
Patel  (1972), 12  are  as  follows: 

THREE-DIMENSIONAL  WORKING  EQUATIONS 

Continuity: 


3u  3v  3w  . 
— — — —  +  — —  +  — —  —  0 
hx3x  3y  3z 


Momentum: 


u  3u 


1  3  /  P 


The  scale  factors  hx  =  hx(y),  hy  =  1  and  hz  =  1  were  used  in  these  equations. 

According  to  Eq.  7.  the  pressure  distribution  is  a  function  of  x  and  z  only,  and 
the  distribution  may  be  determined  from  conditions  at  the  edge  of  the  boundary 
layer.  Thus,  Eqs.  6  and  8  yield 
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For  solutions  to  the  boundary-layer  equations,  the  Reynolds  stresses  are  evaluated 
using  the  Boussinesq  eddy  viscosity  hypothesis 


-  v  w  -i 


Substituting  into  Eqs.  5,  6,  and  8,  the  equations  are  nondimensionalized  by  letting 
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where  c  is  the  chord  length  of  the  airfoil. 
The  working  equations  become 
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For  convenience,  superscript  ~  is  dropped  from  these  equations.  In  the  vortex 
generator  flow  field,  U®  is  replaced  by  Uv,  and  coordinates  x,y,z  are  defined 
accordingly  in  this  region  (Fig.  3). 

On  the  plane  of  symmetry,  w  =  0;  thus  Eq.  (14)  can  be  eliminated.  However, 
for  the  solution,  the  number  of  equations  is  not  sufficient  compared  to  the 
unknowns,  and  an  additional  equation  is  derived  by  differentiating  Eq.  (14)  with 
respect  to  z. 
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At  the  edge  of  the  boundary  layer,  Eq.  15  reduces  further  to 
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The  boundary-layer  displacement  and  momentum  thicknesses  in  streamline  directions 
are  defined,  as 


for  us  =  u  cos  ft  +  w  sin  ft,  (p  =  tan  _  1  , 
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The  flow  on  the  wing  influenced  by  all  the  vortex  generators  is  analyzed  in  a 
small  computational  domain  of  the  wing  span.  The  vortex  generators  are  located 
upstream  of  the  normal  flow  separation  line  for  the  clean-surface  airfoil. 

The  computation  of  flow  characteristics  starts  from  the  forward  stagnation 

point.  The  upstream  laminar  flow  is  assumed  to  be  tripped  by  insect  or  rain 

contamination  and  becomes  turbulent  very  early  at  Re0  =  320.  For  the  purpose  of 

assessing  the  extent  of  the  flow  separation  alleviation,  the  separation  region  has  been 

defined  as  that  downstream  of  the  locations  where  C*  =  0. 

*x 

The  laminar,  transition,  and  turbulent  flow  regimes  are  included  in  the 
computations.  The  extent  of  the  transition  region  is  considered  in  the  computations 
by  multiplying  the  eddy  viscosity  by  the  intermittency  factor  ytr  defined  by  Chen  and 
Thyson  (1971)*3  and  by  Cebeci  and  Smith  (1974). 14 
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where  G  is  the  spot  formation  parameter  defined  by 
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and 

c  =  60  +  4.86  Ms192,  0  <  Ms  <  5 


(for  incompressible  flow  c  =  60). 

An  effort  to  find  a  suitable  eddy  viscosity  model  applicable  to  the  expected 
strong  mixing  of  outer  and  inner  flow  regions  due  to  the  streamwise  vortices  was  not 
successful.  Thus,  the  algebraic  turbulence  model  of  Cebeci  and  Smith  (1974)14  and  its 
extension  was  selected  for  its  simplicity  and  modest  computer  resource  demands.  The 
model  is  described  in  Appendix  B. 


BOUNDARY  AND  INITIAL  CONDITIONS 

For  the  computation,  it  is  assumed  that  the  vortex  line  downstream  of  the  vortex 
generator  is  parallel  to  the  U-direction  and  that  the  vortex  does  not  influence  the 
U-velocity. 

The  boundary  conditions  are 


at  y  =  0,  u  =  v  =  w  =  wz  =  0; 
at  y  =  d,  u  =  U,  w  =  W,  wz  =  Wz; 
at  z  =  0  and  on  the  plane  of  symmetry, 
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The  initial  conditions  of  flow  are  determined  at  the  forward  stagnation  region  of 
wing  leading  edge  using  Hiemenz’s  stagnation  flow  solutions  (Schlichting  I960,13 
pp.  79-83),  as 

^  (Chang  and  Patel,  1975). 16 

FINITE  DIFFERENCE  FORMULATION 

The  Crank-Nicolson  finite  difference  method  is  applied  with  the  grid  structures 
(Fig.  5)  to  solve  the  continuity  and  momentum  equations  of  the  three-dimensional 
turbulent  boundary-layer  flow. 

FORMULATION  AND  ATTACHED  FLOW  GRID  STRUCTURE 

Any  variables  V  and  their  derivatives  are  approximately  given  by  the  following 
equations  (Chang  and  Patel,  1975)16  derived  from  Klinksiek  and  Pierce  (1973):17 

v  =j[v(i+l,j,k)  +  V  (i,j,k)  +  V(i+l,j,k-l)  +  V(i,j,k-l)l  (II 
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If  Ax,  Az,  Ay+  and  Ay_  are  formulated  by 
Ax  =  x(i  +  l)-x(i),  Az  =  z(k)-z(k-l) 

Ay+  =  yO  +  l)-yO).  Ay_  =  y(J)-y(J-l), 

then  the  momentum  Eqs.  13  and  14  can  be  arranged  as  a  tri-diagonal  matrix. 

AO)  u(i  +  l,j-  l,k)  +  BIO)  u(i+  l,j,k)  +  CO)  u(i+  l,j  +  l.k)  =  DIO)  (23) 

AO)  w(i  + 1  ,j  -  1  ,k)  +  B20)  w(i  +  1  ,j,k)  +  CO)  w(i  + 1  ,j  +  1  ,k)  =  D20)  (24) 

Coefficients  AO),  CO),  BIO),  B20),  B30)  are  described  in  Appendix  C. 

On  the  symmetric  plane,  from  Eqs.  13  and  13, 

AO)  u(i  + 1 , j  -  1  ,k)  +  B30)  u(i  +  1  ,j,k)  +  CO)  u(i  +  l,j  + 1  ,k)  =  D30)  (25) 

AO)  wz(i  + 1  ,j  -  1  ,k)  +  B40)  wz(i  +  1  ,j,k)  +  CO)  wz(i  +  1  ,j  +  1  ,k)  =  D40).  (26) 

These  tri-diagonal  matrices  are  solved  by  the  Thomas  algorithm. 

The  velocity  component  v  in  the  y-direction  is  determined  using  the  continuity 
equation  and  applying  the  backward-difference  formula  given  by 

v(i+  l,j,k)  *  v(i+l,j-l,k) 
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Ayj  =  y(j)-y(j-i) 


j  =  2,  3  ...  n  +  1. 


This  equation  is  used  for  the  grid  shown  in  Fig.  S  for  the  solution  of  attached  flow. 


GRID  STRUCTURE  FOR  SEPARATED  FLOW 


For  the  separated  flow,  as  shown  in  the  grid  structure  (Fig.  6),  terms  for  the 
(i+  l,j,k- 1)  points  in  the  region  of  the  separated  flow  region  are  excluded.  Thus,  for 
the  momentum  equation  solution, 


V  =  ^(V(i+l,j,k)  +  V(i,j,k)J 
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and  for  the  continuity  equation  solution, 
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HZ  =  2Az  ~  w  (i.j.k  —  1)  +  w(i.j-l.k)  -  W(i,j-l,k-l)l.  (34) 


The  local  grid  spacings  used  across  the  boundary-layer  thickness  are 
a.  For  laminar  flow 


yj  +  l  =  ^max  ^  +  0  Cm)  j  j,  j-1.  2,  ...,  n. 


Because  the  velocity  changes  rapidly  close  to  the  wall,  grid  size  is  made  small  by 
taking  Cm  <  1.  Boundary-layer  thickness  is  evaluated  by  the  minimum  grid  height 
satisfying 

[ 

For  Cm  =  1,  the  grid  spacing  is  uniform, 
b.  For  turbulent  flow. 

Following  Cebeci  and  Smith  (1974), 14  the  grid  spacing  in  the  y  direction  is 
specified  using  a  geometric  progression  with  the  ratio  of  any  two  adjacent  intervals 


being  a  constant.  Thus, 

Ayj  +  i  =  yAyj 

(37) 

yj  +  l  =  Ay2yy_1  y>l, 

(38) 
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In  y 
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where  Ay2  =  y2-yi-  This  technique  allows  many  grid  points  close  to  the  wall  and 
holds  the  total  points  across  the  boundary  layer  to  a  reasonable  number.  Since  the 
grid  in  the  area  of  the  surface  must  be  smaller  than  in  the  viscous  sublayer,  the 
initial  spacing  was  taken  by  satisfying  Ay2u*/v  <  4.  The  turbulent  boundary-layer 
thickness  was  evaluated  as  in  the  case  of  laminar  flow,  using  Eq.  36.  The  total 
number  of  points  in  the  y  direction  was  in  the  range  of  n  =  24  to  35.  The  step  size 
in  the  x-direction  is  proportional  to  the  boundary-layer  thickness  (typically  Ax  ~  <J). 
The  step  size  in  the  z-direction  is  held  constant  (typically  A z  ~  0.025D). 


RESULTS  AND  DISCUSSIONS 

Numerical  computations  are  carried  out  for  the  vortex  generator  turbulent  flow 
on  a  thick  airfoil  section.  The  airfoil  and  velocity  distributions  were  generated  using 
Schopper’s  (1983)18  airfoil  design  computer  program  which  is  based  on  Lighthill’s 
(1945)19  inverse  method.  The  airfoil  and  velocity  distribution  are  shown  in  Fig.  7. 

Since  such  a  thick  wing  may  be  useful  for  a  high  lift-to-drag  ratio,  low-speed 
cargo  plane,  the  numerical  evaluations  are  performed  for  a  cruising  speed  of  130  ft/s 
and  a  10-ft  chord  length  (Rec  =  8x  107  and  su  =  12.5  ft).  The  angle  of  attack  is  8 
deg  (lift  coefficient  =  1 .03).  The  distributions  of  boundary-layer  shape  factor  H  and 
skin  friction  coeficient  Cf  are  shown  in  Fig.  8  for  a  smooth  airfoil.  As  indicated  by 
the  rapid  increase  of  H  as  well  as  by  the  rapid  decrease  of  Cf  approaching  the 
separation  point,  the  flow  without  vortex  generators  separates  at  x/su  =  0.82,  where 
su  is  the  upper  surface  length  of  the  airfoil.  The  separation  point  is  located  at  the 
position  where  Cf  =  0  and  H  =  2.36.  The  separation  region  is  illustrated  in  Figs.  9, 
10,  and  11. 
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<  0.0005. 


Figure  9  presents  an  overview  of  the  y  =  0  plane  showing  the  vortex  generators, 
the  separation  region,  and  the  computational  domain  of  the  vortex  generator  flow. 

As  shown  in  Fig.  9,  the  separation  line  downstream  of  the  vortex  generator  is  curved. 
Such  non-straight-line  behavior  of  the  separation  line  was  also  observed  by 
Schubauer  and  Spangenberg  (I960).1 

The  extent  of  the  flow  separation  alleviation  due  to  the  vortex  generator  can  be 
assessed  by  comparing  the  separation  line  of  the  vortex  generator  flow  shown  in  Figs. 
10  and  11  with  the  flow  separation  line  without  the  vortex  generators.  The  vortex 
generator  flow  separation  line  is  affected  by  the  position,  size,  and  arrangement  of 
the  vortex  generators,  i.e.,  by  xv,  D,  d,  cv,  h,  ov,  etc.  Since  these  factors  are 
numerous  and  flow  interaction  may  take  place,  it  is  not  an  easy  task  to  interpret  the 
flow  separation  behavior.  Nevertheless,  an  attempt  was  made  to  illustrate  the 
separation  behavior  by  choosing  parameters  D/h  and  D/d  and  by  determining  the 
geometrically  feasible  values  of  cv/h  so  that  no  blade  crossing  would  occur. 
Determination  of  such  cv/h  values  is  described  in  Appendix  D. 

Examples  of  the  separation  lines  are  shown  in  Figs.  10  through  13  for  the 
following  range  of  parameters:  av  =  15  deg,  h  =  0.13  ft  and  0.5  ft,  D/h  =  2-3  and 
10,  D/d  =  2-4,  and  Cy/h  =  1.25  to  8. 

One  example  of  the  investigations  of  effectiveness  of  separation  alleviation  as  a 
function  of  vortex  generator  location  is  shown  in  Fig.  10.  For  these  computations  h 
was  fixed  at  0.15  ft,  and  the  values  of  D/h  and  Cy/h  were  taken  as  10  and  1.25, 
respectively.  Since  good  generator  performance  can  be  obtained  with  a  blade  height 
of  hv  ~  1 .26,  computations  of  6  along  su  indicate  that,  for  generator  locations  of 
xv/su  =  0.4  through  0.7,  the  fixed  generator  height  value  may  serve  to  allow 
assessment  of  the  effect  of  generator  location  on  generator  performance.  As  shown 
in  Fig.  10,  the  separation  line  is  located  farther  downstream  as  the  generator  is 
moved  upstream.  The  figure  also  shows  that  performance  improves  as  the  spanwise 
spacing.  D/d,  increases. 

Figure  1 1  presents  separation  line  results  for  the  generator  fixed  at  xv/su  =  0.4 
and  various  values  of  the  parameters  D/h,  D/d,  and  cv/h.  For  this  study  h  =  0.5  ft 
is  assumed,  and  although  this  value  is  larger  than  1 .2d  and  the  results  are  strongly 
affected  by  turbulence  mixing,  the  value  appears  to  be  suitable  to  evaluate  the 
correlation  of  geometric  and  position  parameters.  Figure  11  also  shows  that  the 
separation  point  moves  downstream  along  the  symmetry  axis  A  and  moves  upstream 
along  the  symmetry  axis  B  ahead  of  the  separation  line  which  occurs  without  the 
vortex  generator. 

Wang  (1970)20  found  that,  for  the  three-dimensional  boundary  layer  near  the 
plane  of  symmetry  of  a  spheroid  at  an  angle  of  attack,  at  least  near  the  symmetry 
plane,  separation  occurs  upstream  of  the  streamwise  separation  point  defined  by  the 
location  at  which  the  surface  streamlines  from  each  side  converge  and  meet 
tangentially.  Wang  conjectures  that  this  phenomenon  is  attributable  to  the  cross-flow 
separation  (i.e.,  due  to  the  reversal  of  v  velocity).  Further  investigations  of  cross-flow 
separation  of  the  vortex  generator  may  clarify  this  phenomenon. 

Since  the  extent  of  the  upstream  region  of  separation  beyond  the  separation  line 
without  the  vortex  generator  is  small,  the  vortex  generator  is  effective  with  larger 
cv/h  values,  although  the  effectiveness  does  not  increase  much  and  remains  nearly 
unchanged  with  values  of  cv/h  >  6.  It  should  also  be  noted  that,  with  equal  values  of 
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D/h  and  Cy/h,  the  separation  region  is  reduced  as  D/d  increases,  indicating  higher 
efficiency  of  separation  alleviation  with  larger  D/d. 

Figure  12  presents  shape  factor  distributions  as  a  function  of  surface  distance 
and  parameters  D/d,  D/h,  and  z  when  the  vortex  generator  is  fixed  at  xv/su  =  0.4. 
Over  most  of  the  acceleration  portion  of  the  airfoil,  x/su  <  0.55,  the  shape  factor  is 
nearly  constant  at  H  ~  1.5.  In  the  airfoil  deceleration  region,  the  shape  factors  show 
earlier  and  steeper  growth  with  increasing  z,  that  is,  with  increasing  distance  from  the 
symmetry  axis  A.  Along  the  symmetry  axis,  the  results  are  nearly  independent  of 
D/d,  D/h,  and  z. 

Figure  13  shows  streamwise  skin  friction  distribution  when  the  generator  is 
located  at  xv/su  =  0.4.  Here  Cf  decreases  rapidly  in  the  deceleration  region  and 
reaches  the  value  of  zero  at  the  separation  point.  With  increase  of  z,  separation 
occurs  earlier.  As  seen  in  Fig.  13c,  for  D/h  =  3  and  4  separation  is  completely 
suppressed  up  to  the  trailing  edge. 

Additional  study  of  the  computed  results  (not  reported  herein)  may  establish 
optimum  conditions  of  the  vortex  generator  flow  to  alleviate  the  flow  separation.  It 
would  be  of  interest  to  investigate  the  extent  of  separation  alleviation  by  installing  a 
staggered  double  row  of  vortex  generators.  The  analysis  could  be  improved  by 
applying  the  newest  turbulence  model  of  Cebeci  et  al.  (1987)21  for  two-dimensional 
and  three-dimensional  turbulent  flow  involving  flow  separation.  Since  the  vortex 
generators  cause  drag,  a  drag  analysis  of  the  vortex  generator  also  would  be  of 
interest. 


CONCLUSIONS 


A  numerical  turbulent  flow  analysis  of  counter-rotating  vortex  generators 
attached  to  a  large  aspect  ratio  wing  has  been  performed  for  the  case  of  Rec  =  8  x 
107.  The  Cebeci  et  al.  turbulence  model  has  been  used  under  the  assumption  that 
transition  occurs  near  the  leading  edge  at  Ree  =  320.  The  following  conclusions  have 
been  reached: 

1.  With  an  increase  of  Cy/h,  the  separation  point  moves  downstream;  along  the 
symmetric  axis  of  the  vortex  generator  (z  =  0).  Along  the  symmetric  axis 
between  the  vortex  generators  (z  =  D/2),  the  separation  point  is  located  a 
short  distance  upstream  of  the  separation  line  without  the  vortex  generator. 
Nevertheless,  the  effectiveness  of  the  vortex  generator  in  alleviating  flow 
separation  is  evident.  When  cv/h  is  large  (greater  than  about  4),  separation  is 
completely  suppressed  in  the  axis  region  of  the  vortex  generator.  Beyond  a 
certain  limit,  perhaps  cv/h  >  6,  the  effectiveness  of  the  vortex  generator  does 
not  increase  appreciably. 

2.  Separation  is  more  effectively  suppressed  as  the  vortex  generator  is  moved 
upstream  and  as  the  spacing  between  the  pairs  of  counter-rotating  vortices  is 
increased. 

3.  The  behavior  of  the  distributions  of  H  and  Cf  along  the  vortex  generator 
symmetric  axis  is  nearly  identical  for  all  cases  of  separated  flow  computed 
and  is  not  affected  by  D/d,  D/h,  and  cv/h,  but  depends  on  the  spanwise 
distance  z.  With  a  decrease  in  z,  separation  occurs  later,  and  in  some  cases 
separation  is  prevented. 


13 


j. 


*  I 


l\\  I  \ 


I 

/ 


\ 


SYMMETRIC  AXIS  B 


FLOW  DIRECTION 

I  * 

'I' 

l\ 

SYMMETRIC  AXIS  A 

SYMMETRIC  AXIS  A 


/  \ 


B 

y 

*  A 

rv 

\ 

VORTEX  STRENGTH:  T 

r  ^ 

L _ _ _ ^ 

XVnSVsNXVvSxxxVnnnVnv 

c 

^WWWWW  \\\\\\\\\vCx — 

D 

»  z 


Fig.  3.  Vortex  generator  arrangement. 


Fig.  4.  Coordinate  system. 
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Fig.  7a.  Airfoil  configuration. 
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Fig.  7b.  Velocity  distribution  on  airfoil  surface  (o  =  8°). 

Fig.  7.  Sketch  of  configuration  and  velocity  distribution  of  airfoil. 


Fig.  8a.  Shape  factor  H  along  main  stream  direction. 


Fig.  8b.  Skin  friction  coefficient  Cf  distribution. 

Fig.  8.  H  and  Cf  behavior  in  absence  of  vortex  generator  (o =8°). 


Fig.  9a.  Separation  line  for  weak  alleviation. 
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Fig.  10c.  D/d  =  2. 
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Fig.  10.  Separation  line  in  parameter  of  xv/su. 
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Fig.  12*.  D/d  =  2,  cT/b  =  2. 
hape  factor  H  curve  in  parameter  of  z. 


Fig.  12.  (Continued) 
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APPENDIX  A 


APPROXIMATE  EVALUATION  OF  VORTEX  STRENGTH  OF 
VORTEX  GENERATOR 

The  two-dimensional  approximate  analysis  of  Jones  (1957)3  is  used  to  evaluate 
the  vortex  strength  of  the  vortex  generator;  more  details  are  given  in  Ref.  3.  The 
induced  angle  of  attack  to  the  vortex  generator  blade  is  W/Uv.  Thus  the  effective 
angle  of  attack  to  the  blade  at  any  point  P  is 


♦v  av 


w 

U„ 


The  lift  at  any  span-wise  element  dz  is 

W„ 


d  Lv  =  a  i  *v  - 


icv-Uv2  dz 


U„  ( 


where  a  is  the  lift  slope. 

For  convenience  the  mean  value  Wm  is  taken  in  order  to  apply  the  simple  horseshoe 
vortex  system  of  lifting  line  theory. 

The  Kutta-Joukowsky  theorem  gives 


<3  Lv  =  q  uv  T  dz 
Hence, 

W 


r  = 


Letting 


m 


Uv 


}cvUv. 


where 


DUV 


cv/D 


a  s 

2  i  “v 


uriCi’ 


and  D  is  the  distance  between  the  vortex  generators  shown  in  Fig.  3,  and 
W„ 


[A 
D  ’ 


the  previous  equation  becomes 


y  = 


--  +  a 


aci 
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The  total  strength  of  the  trailing  vortices  behind  any  generator  is  assumed  equal 
to  the  total  circulation  around  that  generator,  i.e.,  around  blade  of  height  h. 
Therefore,  if  is  used  instead  of  D  and  y'  is  defined  by 


r  =  y'  cv  Uv  =  Y  DUV. 


D  1 

y  =  r  7-  =  /  7- 


y  = 


r  = 


aci 

Os, 


+  A 


—  +  A  C‘ 
ac, 

,  Wm  D  Wm  D 

A  =  “f—  =  -7~Z  '  Cy 


Cl 


—  -c  -U  =  _gv  cv  Uv 


a  +  C>A 


=  fv  _J_  /*h  _W 

‘  Cj  h  j  r 


dy 


for  a  given  x  and  z.  Hence, 


r  = 


ffv-cv  Uv 


2  cv  /*  W  J 

7  +  rdy 


(a  ~  2n) 


Therefore 

r  = 


drv  cv  Uv 


1  cv  /*  W  J 

^+hV  T dy 


or 


r  ~ 


av  h  Uv  ARV 


-  +  AR, 


I 


where 


APPENDIX  B 

CEBECI  AND  SMITH  TURBULENCE  MODEL 
Cebeci  (Refs.  14,  22-24)  proposed  the  following  eddy-viscosity  model: 

£  =  (>cy)2  [l  -  exp  (~  X^)]2/(|y)2  +  (|y)2  ;  0  <  y  <  yc  (Cebeci  (1975]22) 

£  -  a  |  f  (Us  -  Vu2  +  w*)  dy  |  [l  +  5.5  ^  6j~ 1;  yc  <  y.  (Cebeci  [1970]23) 

Since  t  is  continuous,  yc  is  determined  from  these  two  equations,  and  the  upper  limit 
satisfies  the  law  of  wall  y*  =  y  u*/v.  A*  is  the  Van  Driest  dampling  factor  given  by 


A*  =  Aq/N 


where 


w  i/T - iTa~^«  3p  v  (Cebeci,  Smith  [1974],14 

N  =  VI  -  11.8  p*  p  =  -  _  pp  2J4>  2]6  21?) 

o  _ 

is  the  pressure  gradient  along  the  stream  line 

r^- 

u*  =  /  —  =  wall  shear  velocity 
”  Q 

x  =  0.4  +  — — - ,  o  =  Ree10~3  >  0.3  (Cebeci,  Smith  [1974], 14  p.  222). 

1  +0.49  o2 


Aq  =  26  + 


q  =  0.0168 


1  +o2 


Ree  >  6000 


0.0168  x  1.55 


(Cebeci  [1973], 24 


425  *Kee  <  6000  Cebeci>  Smith  [I974]f14  p.  224). 


a  =  10“ 3  (194.8- 128.6 x log  Refl+ 30.925  x  log  Re02 


-  2.475  x  log  Re03),  Re0  <  425 


(Chang,  Patel  [1975]  *6) 


where 


Rea  * 


'‘WWW 


APPENDIX  C 

COEFFICIENTS  OF  MOMENTUM  EQUATIONS 

Ax  =  x(i+l)  -  x(i),  Az  =  z(k)  -  z(k-l) 

Ay^  =  y(j+  1)  -  y(j).  Ay_  =  y(j)  -  yO~l) 
hx  Ax  hx  Ax 

,  A _ _ _  _l  _ _ _ _ — _ 

+  =  Ay+  (Ay+  +Ay_)  ’  '  _  Ay_  (Ay+ +Ay.) 

v  d  +  Ay  _  V  d  _  Ay  +  hxAxWz 

r'  =  20  ’  f2  =  20  ’  =  20 

Sj  =  [1  +t(i  +  l,j  +  1/2, k)]  d  +  /u 

s2  =  [1  +  e(i,j -  1/2, k)]  d_/u 

s3  =  [1  +c(i,j  +  l/2,k)]  d  +  /u 

s4  =  [1  +e(i,j  -  1/2, k)]  d_/u 

s5  =  [1  +c(i  +  l,j  +  1/2, k-  1)]  d  +  /u 

s6  =  [1  +£(i  +  l,j  -  1/2, k-  I)J  d_/u 

S7  =  [1  +t(i,j  +  1/2, k  —  1)]  d  +  /u 

s8  =  [1  +  c(i,j  -  I/2,k-  1)]  d_/0 

W  hx  Ax 
1  ~  TT  Az 

AO)  =  -  r2  -  s2,  CO)  =  rj  -  s, 

B]0)  =  1  +  r2  -  rj  +  s2  +  Sj  +  t 

B20)  =  1  +  r2  -  r,  +  s2  +  s,  +  t 

B30)  =  1  +  r2  -  r,  +  s2  +  s, 

B40)  =  1  +  r2  -  r,  +  s2  +  sj  +  r3 

D,0)  =  4r[U(i  +  l,k)  -  U(i,k)  +  U(i+l,k-l)  -  U(i,k-1)] 

_  HxA^W  [U(i +  l,k)  -  U(i,k)  -  U(i,k-  1)  -  U(i+l,k-l)] 


+  (r2  +  s4)  u(i,j-  l,k)  +  (1 +  r, -r2-s3-s4  +  t)  u(i,j,k) 

-  (-ri  +s3)  u(i,j  +  l,k)  +  (rj  +  Sft)  u(i,j- l,k- 1) 

+  (1  +r,  -r2-s7-sg-t)  u(i,j,k  -  1)  +  ( -  r,  +  s7)  u(i,j  +  l,k  -  1) 


D2(j)  = -=r[w(i  +  l,k)  -  w(i,k)  +  w(i+l,k-l)  -  w(i,k-l)J 


hx  Ax  W 


(w(i+l,k)  +  w(i,k)  -  w(i  +  l,k-l)  -  w(i,k-l)] 


-  (r2-s4)  w(i,j  -  l.k)  -  (1 +r, -r2-s3-s4-t)  w(i,j,k) 

+  (-r,  +  s3)  w(Lj  -  l,k)  +  (r2  +  St)  w(i,j-l,k-l) 

+  (1  +  r,  —  r2  —  s7  —  s8  —  t)  w(i,j,k -  1)  +  (s7-r,)  w(i,j  +  l.k-  1) 

D3G)  =  -=-[U(i  +  l,k)  -  U(i,k)]  +  (r2  +  s4)  u(i,j-  l.k) 


+  (1 +r,-r2-s3-s4)  u(i,j,k)  +  (-r,  +  s3)  u(i,j  +  l,k) 


+  (r2  +  s4)  w2(i,j  -  l.k)  +  (1 +r,-r2-s3-s4-r3)  wz(i,j,k) 
+  (-r,  +  s3)  wz(i,j-l,k) 


Superscript  —  is  space  averaged;  u,  v,  w,  and  wz  are  assumed  to  be  known  by 
projecting  the  velocity  profile  at  the  station  (i,j,k)  forward  to  station  (i  +  l,j,k)  for  all 
y  and  k  according  to 


V(i-»-  l,j,k) 


V(i  +  l,n  +  l,k) 
V(i,n  +  l.k) 


V(i,j,k). 


APPENDIX  D 

GEOMETRICAL  ARRANGEMENT  OF  VORTEX  GENERATORS 

Schopper*  investigated  the  problem  of  determining  the  maximum  chord  length  cv 
which,  for  a  given  D/h  and  D/d,  does  not  result  in  the  blades  crossing.  The  term 
cv/h  is  evaluated  as  a  function  of  (D/h)/(D/d)  =  d/h. 


a.  Self-crossing 


For  a  given  d/h,  self-crossing 
occurs  when 


(cv/h). 


d/h 


max  sin  a 

For  no  self  crossing, 

cv/h  <  (cv/h) 
v  v  max 


The  condition  for  no 
interference  between  pairs  is 

D'  dj  +Cy  sin  a 

T> - K - 


1  d  -  cv  sin  a 


d» 

T  "  2 


D-dj  -cv  sin  o 
h  h 


Then 


(40) 


D-dt  -cv  sin  a  dj  +  cv  sin  a 


or 


Cy  .  D-d  1 

h  h  sin  a 


or 


D  ^  d  +  cv  sin  a 
h  h 


•Private  communication  with  M.R.  Schopper,  David  Taylor  Research  Center. 


For  no  interference, 


£v  <!  1 

h  h  sin  a 


D 

d 


-r  -  1 


(41) 


From  Eqs.  40  and  41  and  D/d  >  2,  cv/h  is  determined  by  Eq.  40.  Examples  of 
(cv/h)max  are  shown  in  the  following  table  for  several  values  of  D/d  and  D/h  at  av 
=  15  degrees,  i.e.,  1/sin  av  =  3.864: 


D/d 

D/h  =  2 

D/h  x  4 

D/h  =  5  j 

D/h  =  10 

d/h 

(cv/h>max 

d/h 

^Cv^max 

d/h 

^v^^max 

d/h 

(c./h^x 

d/h 

<Cv/h>max 

2 

1 

3.864 

1.5 

5.796 

2 

7.728 

2.5 

9.66 

5 

19.32 

3 

2/3 

2.576 

1 

3.864 

E2M 

5.15 

9EKUH1 

IBM 

”15.06 

4 

1/2 

1.932 

3/4 

2.898 

1 

1004 

"“5.66  ' 

5 

2/5 

1.546 

3/5 

2.318 

4/5 

3.09 

1 

i  3.864 

rrn 

7.728 
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